Stationary localized modes in the quintic nonlinear Schrodinger equation with a 
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We consider the localized modes (bright solitons) described by one-dimensional quintic nonlinear 
Schrodinger equation with a periodic potential. In the case of attractive nonlinearity we deduce 
sufficient conditions for collapse. We show that there exist spatially localized modes with arbitrarily 
large number of particles. We study such solutions in the semi-infinite gap (attractive case) and 
in the first gap (attractive and repulsive cases), and show that a nonzero minimum value of the 
number of particles is necessary for a localized mode to be created. In the limit of large negative 
frequencies (attractive case) we observe "quantization" of the number of particles of the stationary 
modes. Such solutions can be interpreted as coupled "Townes" solitons and appear to be stable. 
The modes in the first gap have numbers of particles infinitely growing with frequencies approaching 
one of the gap edges, which is explained by the power decay of the modes. Stability of the localized 
modes is discussed. 
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I. INTRODUCTION 

The quintic nonlinear Schrodinger (QNLS) equation 

iu t + u xx — V{x)u — <j\u\ a u = , (1) 

where either a = +1 (the repulsive case) or a = —1 
(the attractive case), and V(x) is an external potential, 
is a standard model of the nonlinear physics. This model 
can be viewed as the evolution equation u t — iSH[u]/5u* 
generated by the Hamiltonian 



H[u] = f 

J —t 



[\u x \ 2 + V(x)\u\ 2 + f(\u\ 2 )]dx (2) 



where the nonlinearity /(|u| 2 ) has the second order term 
of the Taylor expansion equal to zero, i.e. /(|m| 2 ) = 
(a/3)\u\ 6 + 0(|w| 8 ) (here we took into account that 
/(0) = 0, necessary for convergence of the integral, and 
dropped the unessential term /'(0)|u| 2 which results in 
adding a constant to the Hamiltonian). In particular, 
this is the case of a slow envelope of a small amplitude 
wave in a medium with weak dispersion and nonlinearity 
having specific dependence on the wave intensity. 

Recently, renewed interest in the model Q was orig- 
inated by the suggestion of its use for description of a 
one-dimensional gas of impenetrable bosons In the 
later studies Q , limitations of the mean-field approach 
have been discussed. More specifically, it has been shown 
that the model fails at relatively small number of atoms 
and in description of the gas dynamics. Nevertheless, 
it has been rigorously proven in that in the limit of 
large number of particles, there exist a domain where the 
stationary version of Eq. [see 10} below] , correspond- 
ing to the energy functional given by the integral @ 
with f{\u\ 2 ) = (7r 2 /3)|w| 6 (provided V(x) has a homoge- 
neously growing part), does describe the ground state of 
the Tonks- Girardeau gas Q. 



Let us now assume, without loss of generality, that 
V(x) is a 7r-periodic function, V(x) = V(x + n). This 
means that the energy, and in particular V(x), is mea- 
sured in units of the recoil energy Er = H 2 tt 2 / (2md 2 ), 
where m is a mass of an atom and d is the lattice con- 
stant, and time and spatial coordinate in the dimension- 
less Eq. Q are measured in the units Er/H and d/ir, 
respectively. Then, we observe that the number of real 
atoms of the Tonks-Girardeau gas, which we denote as 
AT, is connected with the integral of motion of the di- 
mensionless model Q N = J \u(x,t)\ 2 dx (hereafter we 
drop limits in integrals over the real axis) by the formula 
TV = N/ir (for the sake of brevity below N is also referred 
to as a number of particles). 

Next, we recall that in the case of the cubic nonlinear 
Schrodinger equation, where /(|it| 2 ) oc |u| 4 , with a peri- 
odic potential having amplitude of the unity order, the 
quantity N for a gap soliton is of order of a few units [f| , 
and that localized solutions of the QNLS equation re- 
ported in earlier publications @ have also N of order of 
one (see also the results presented below) . The respective 
modes do not belong to the domain of the applicability of 
the model |QJ in the theory of the Tonks-Girardeau gas 
(i.e. do not satisfy the necessary condition N 3> 1) |15| . 
This rises a natural question: can localized modes (or gap 
solitons) of the model Q have any physical relevance to 
the theory of the Tonks-Girardeau gas? 

In the present paper we will give a positive answer 
to this question, which will be based on the argument, 
that the QNLS equation with a periodic potential allows 
storage of an infinite number of particles, i. e. allows 
lcalized mode solutions corresponding to the limit N 3> 
1. 

Meantime, the model appears also in the quasi- 
one-dimensional limit of the generalized Gross-Pitaevskii 
equation when two-body interactions are negligibly small 
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(what can be achieved, say by means of the Feshbach 
resonance) and the chemical potential is approached by 
/i w giv? with the coefficient 32 depending on the three- 
body interactions and n being the atomic density @, 
(see also |8|). In that case, the total physical number of 

particles is computed as TV = \J 2^m\^ r\^ ( nere a ± m a 
transverse linear oscillator length, which must be much 
bigger than the scattering length a s : a± ^> \a s \). Also in 
this case a =sign (32)- 

The homogeneous case V(x) = of Eq. has been 



extensively studied (see e.g.|9|)- It is known that in the 
attractive case Eq. Q describes two specific scenarios of 
the evolution. The first scenario corresponds to collapse 
when a smooth solution ceases to exist and its amplitude 
infinitely grows in a finite time. The second scenario 
is dispersion of a solution, when amplitude of initially 
localized solution decays and the area of its localization 
grows. These two regimes are "separated" by the soliton 
solution (see also Sec. II below) 



u(t, x) 



3 l/4(. 



a/4 



y/ cosh(2\/— lux) 



(3) 



Here a; is a dimensionless real constant parameterizing 
the solution. As it is clear solution J3J exists for u> < 0. 
Bellow lj will be referred to as frequency, while in the con- 
text of the background solution within the framework of 
the mean-field theory of the quantum gases 0, 0] w de- 
termines the chemical potential fi, given by fi = Erco/H. 

In the repulsive case neither bright soliton solutions 
nor collapse regime can occur. 

Presence of the periodic potential V(x) in the model 
greatly enriches its behavior. In Ref. |6j it has been re- 
ported the existence of unstable solutions which allow 
one to characterize the delocalizing transition ^(| and 
respectively the existence of the minimal and maximal 
values of the number of particles corresponding to exis- 
tence of localized solutions beyond which the solution is 
either collapsing or dispersing. 

It turns out that the earlier results can be signifi- 
cantly generalized. In particular, in the present paper 
we show that in the attractive case and for frequencies 
laying in a semi-infinite gap solitary waves represent an 
infinite number of branches of the solutions which can be 
parametrized by the number of particles. We also show 
that the presence of the lattice significantly affects prop- 
erties of the model, on the one hand changing sufficient 
conditions for the collapse and on the other allowing one 
to store an arbitrary number of particles either in a form 
of a soliton chain or in a single soliton. We compute the 
bifurcation diagram for the lowest branches of such solu- 
tions and describe a phenomenon of " quantization" of the 
number of particles in the limit of large negative frequen- 
cies in the attractive case. We also describe families of 
the gap-soliton solutions showing great differences with 
the case of the cubic nonlinear Schrodinger equation. 



II. STATEMENT OF THE PROBLEM 

We assume that V(x) is bounded, |V(x)| < Vo (Vo is a 
positive constant), 7r-periodic, and even, V(x) = V(—x), 
with a local minimum placed at x — 0. We seek 
for stationary localized solutions in the form u(x, t) = 
exp(— iu>t)<f>(x) where w is a frequency. The function 4>{x) 
can be considered real and satisfying the equation 



<\>xx ~ <T<P 5 + [(J- V(X)] = 



(4) 



with zero boundary conditions: lim^i^oo 4>(x) — 0. A 
localized solution decays at infinity, thus reaching spa- 
tial domains where the nonlinearity becomes negligible. 
Therefore its asymptotics are determined by the linear 
eigenvalue problem 



d 2 <p aq 

dx 2 



+ [tO a (q) - V(x)\ (faq = 



(5) 



-) ,.(+) 



which has a band spectrum [lfj uj a (q) € 

with a — 1,2,... referring to a band number and "+" 
standing respectively for the upper and lower 



and 



boundaries of the a-th band. ip aq are Bloch functions 
and q is the wave number in the reduced Brillouin zone, 

q € [—1,1]- An interval (cj^~\ w^jj represents the a's 

gap. The lowest gap, which we designate by a = 0, 

is semi-infinite, ^— 00, lu[ ^. Since the asymptotics of 

localized solution of Eq.QJ is determined by (0, there 
exist no localized solutions of Eq. Q for to lying in a band 
and only for to lying in a gap the possible. 




FIG. 1: The structure of the gaps for Eq.@, V(x) = 
— Vbcos2a;. The triangular area in the lower part of the fig- 
ure is the area ui < Vo where no localized solutions exists in 
repulsive case (see the text). 

As we already mentioned, an important characteristics 
of the solution is the integral N, which in low-dimensional 
condensate applications is proportional to the number of 
particles of the condensate, therefore below we refer to 
N as to a number of particles. 

In the attractive case and in the absence of the pe- 
riodic potential, V(x) = 0, the number A* = \/3ir/2 
separates the solutions with N > A» which may col- 
lapse and solutions with N < N* for which no collapse 
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can occur . Specifically, TV, corresponds to the exact 
solution @ and the frequency, uj < 0, can acquire an 
arbitrary negative value. This situation is similar to one 
discovered in the context of the nonlinear optics for the 
2D cubic Schrodinger equation having the stationary so- 
lution, separating collapsing and dispersing pulses. That 
solution was discovered in 1964 and today is referred 
to as Townes soliton. Using this analogy, the stationary 



solutions of with 



uj < uj 



(-) 



, i.e. in the semi-infinite 



gap, having the smallest number of particles at given uj, 
will be referred to as Townes solitons (TSs). Thus, by 
TSs we call the one-parametric family of solutions cor- 
responding to the lowest branch in the semi-infinite gap 
(see the discussion below and Figs. |21 and |3Jl. In the limit 
V(x) = TSs are transformed in the solitons given by 
(0. If V(x) ^ stationary localized solutions of 

with uj lying in a finite gap, uj G ^a+ij > mav ex ~ 

ist. Also, the new solutions, with the number of particles 
bigger than the number of particles of TSs, may exist 
in the semi-infinite gap. All such solution disappear in 
the limit V{x) = 0, and therefore we will call them gap 
solitons. 



III. SUFFICIENT CONDITION FOR COLLAPSE 
IN THE ATTRACTIVE CASE (a = -1) 

From simple intuitive arguments it follows that the sta- 
bility of an initially strongly localized wave packet, (i.e. 
localized on distances much smaller than the lattice pe- 
riod), should not be significantly affected by the lattice. 
Indeed, boundness of J V<j> 2 dx and use of the known ar- 
guments for obtaining iV* 0, allow one to conclude 
that the presence of a periodic potential does not affect 
the fact that no collapse can occur if N < . 

To deduce a sufficient condition for the collapse for 
Eq.Q we introduce two momenta, 



y(t) 



x 2 \u(x, t)\ 2 dx , 



z{t) = Im J xu(x,t)u x (x,t)dx 



and the energy functional 



E=\j (\u x \ 2 -\\u\ G + V\u\^dx. 



(6) 



which is a first integral of Eq.QJ. Next, following the 
standard steps , and assuming that the solution exists 
for < t < T, we obtain 



d 2 y A dz 
d^ 



dr f 

-4— = 16£ - 4 / (2V + xV x )\u\ 2 dx (7) 



Let us show that the conditions 

z(0)>0 and r, = -AE - V iy 1/2 (0)N 1/2 > , (8) 



where E = E + V N/2 and V% = sup.,, |V^(a;)|, are suf- 
ficient conditions for collapse in presence of periodic po- 
tential. 

First, one can show that y(t) is a decreasing positive 
function, i.e. < y(t) < y(0). Indeed, define v{t) = 
y{t) — 2/(0) and obtain from Q 



dv(t) 
dt 



4z(0) + i 



v(s)ds = F(t) 



(9) 



Then dv/dt < for all t £ [0, T). Indeed, let us assume 
that there exists i» < T such that dv(t*)/dt = 0, while 
dv/dt < for t < f*. From the definition of F(t) it follows 
that F(0) < 0, dF/dt = -An + v/2 < for t € [0,t,) 
and, consequently, F(t*) < F(0) < 0. This contradicts 
the condition F(t#) > which follows from (0. Thus, 
v(t) decreases at t £ [0, T). In addition, Eq. J7J) implies 
that 

y(t) - 2/(0) < (8£; + 2V l y 1 / 2 {0)N 1/2 )t 2 ~ 4s(0)* (10) 

Thus we have outlined the proof of the following state- 
ment: 

Theorem - Let conditions be satisfied. Then. 
y{t) < 2/(0) — 2r\t 2 and collapse occurs at finite time 
T<(?X0)/2 ?7 ) 1/2 . 

Comparing the sufficient conditions (JSJ) with the con- 
ditions for collapse for homogeneous case V\ — we con- 
clude that in the presence of a periodic potential the col- 
lapse occurs at smaller values of energy E < 0. This re- 
sult has simple physical explanation. Indeed, for atomic 
distributions involving several local minima, the respec- 
tive maxima act act as a repulsive force diminishing the 
effect of the inter-atomic attraction. 



IV. A MINIMUM NUMBER OF PARTICLES 
FOR A LOCALIZED SOLUTION 

When N < iV*/v3, multiplying @ by <j>, integrating, 
and using the Gagliardo-Nirenber g in equality in the form 
J(j) 6 dx < {3N 2 /N 2 ) J (j) 2 x dx (see p|) we obtain the es- 
timate 



3 N 2 

( j } ldx<(cj + V )N[l + - — (<7-l) 



(11) 



This inequality has the following nontrivial consequence: 
For the stationary solutions of Eq. 0) the limit of a small 
number of particles, N — > 0, if available, implies small- 
ness of the amplitude of the function, sup x \4>\ — > 0. 
The stationary solutions at uj < —Vo (if any) have 
N > Njy/3. 

This statement follows immediately from (|11|) and the 
inequality (sup x </>) 4 < AN J (f> 2 dx. Thus, in order to es- 
tablish the existence of a lower bound for the number of 
particles necessary for creation of a stationary localized 
solution, we address the small amplitude limit of Eq. . 
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The respective solutions bifurcate from the linear 
band, and thus can be described in terms of the effective 
mass approximation p|. Thus we look for a solution of 

||TJ in a form ijj — t x l 2 i\)\ + e 3 / 2 V>2 H f where e< 1 and 

are functions of slow variables x p — e p x and t p — e p t. 
Then, following the standard procedure, we obtain that 
ipi = A(xi, t2)(f a (xo) exp(— iu a to), where u> a stands ei- 
ther for u>a (if a = 1) or for oj^X (if a — —1) and ip a 
stands for the respective Bloch function [see J5J)]. The 
slow envelope A(x\,t2) solves the QNLS equation 



.8A 



1 d 2 A 



dt 2 2M a dx\ 



a Xa \A\ i A = {) 



(12) 



where Xa = Jq \<fia(x)\ e dx is the effective nonlinearity 
and MqT 1 = d 2 uj a (qo)/dq 2 is the effective mass (go = 
and qo = 1 for the center and the boundary of the 
Brillouin zone, respectively). The soliton solution of H12|) 
is known: 



A = 



M 4 



cosh I exy/8\M a 



(13) 



Thus e 2 can be interpreted as a frequency detuning to 
the gap, e 2 = \u> — w a |, and the number of particles 



J \A\ 2 dx = 



2^2\ Xa M a \ 



(14) 



does not depend on the detuning, but is determined 
only by the lattice parameters. This contrasts with 
the situation which takes place in the case of nonlin- 
ear Schrodinger equation with cubic nonlinear term: the 
number of particles i V there te nds to zero together with 
the detuning: N oc y/\u) a — ui\ (see e.g. 0). 

The above arguments prove the existence of the min- 
imal number of particles necessary for the creation of 
stationary localized mode described by Eq.lJIJ. Indeed, 
if N — > would be possible, it necessarily should be 
reached at solutions . The latter, however, does not 
allow such a limit. At the same time, according to the 
theory presented, the obtained estimate for N itself needs 
not be the smallest number of particles of the stationary 
solutions. And it is not, as is seen from examples of nu- 
merical studies shown in Figs. El and [S] where, like in 
all other simulations, we use the potential 



V(x) = -Vocos(2x) 
and show only the three lowest branches. 



(15) 



V. SEMI-INFINITE GAP, TOWNES SOLITONS. 
A. The attractive case. 
Fig. |21 shows the behavior N vs u in the semi-infinite 
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gap in the vicinity of its upper bound u}\ 



(-) 



-0.937. 



FIG. 2: Number of particles vs frequency u> for even TSs at 
Vo = 3 (the central box). In the side boxes explicit forms of 
solutions are shown. The box letters correspond to the points 
on the branches: (A) for u) = -1.2, (B) for w = -1.1, (C) 
for lo = -0.9374 and (D) for u = -0.95. u, fa 1.0005 is the 
bifurcation point of the second and third lowest branches. 



The number of particles corresponding to TSs (the low- 
est branch) achieves its non-zero minimum, showing non- 
monotonic behavior. At large frequency detunings to- 
ward the gap new branches of the solution appear. Two 
of them bifurcate at some point u> = a;*. Fig.[5]also shows 
that increase (decrease) of the number of particles in a 
TS necessarily lead neither to collapse nor to spreading 
out of the solution. Indeed, simultaneous change of the 
frequency and the number of particles, corresponding to 
motion along the respective branch, simply modifies the 
stationary solution, what is expressed in the reshaping of 
the mode: when frequency approaches the gap a mode 
broadens and its amplitude decreases (c.f. boxes B, C 
and D). 

Existence of the local minimum on the lowest curve 
N(ui), at a frequency shifted from the edge toward the 
gap, does not contradict the small amplitude approxima- 
tion which gives (|14|) . thus predicting N to be constant. 
Indeed, the small parameter e 1 / 2 ~ 0.1 corresponds to 
the frequency detuning \u — ui a \ ~ 10 -4 i.e. is extremely 
small, therefore the region of validity of asymptotic ex- 
pansions is also very narrow. This is illustrated by the 
panels B, C and D in Fig. |3 showing that reshaping of 
a soliton (from the envelope soliton in C to an intrinsic 
localized mode in B) occurs already at quite small fre- 
quency detuning. It is worth it to compare this result 
with the one known for the cubic nonlinear Schrodinger 
equation, where the number of particles is a monotonic 
function of the frequency detuning Q: A a cx \u> — LJa] 1 ^ 2 
and the small parameter is s instead of e 1 ! 2 obtained 
above. 

As it is illustrated by the panels B,C, and D, the low- 
est branch describes a family of one-soliton solutions. In 
a similar way the upper branches can be associated with 
multi-soliton solutions. As an example in panel A of 
Fig. |3 we show a two-soliton mode corresponding to the 
second lowest branch. In the point uj = shown in the 
central panel, this mode bifurcates with a three-soliton 
solution (not shown here) corresponding to the next up- 
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per branch. In the same manner one can describe the 
whole set of infinite number of upper branches. 




4 -2 2 4 



FIG. 3: Number of particles vs frequency for even TSs, at 
Vo — 1 (solid line), Vb = 3 (dashed line), and Vo = 5 (dotted 
line), corresponding to first, second and third branches shown 
in Fig|5| In the boxes (a) , (b) , and (c) we show explicit shapes 
of the solutions for Vo = 3 and u = —12.175. 



lies below the critical value AT*. This means that, unlike 
the homogeneous QNLS equation > small enough per- 
turbations of stationary solutions in presence of a peri- 
odic lattice do not result in collapsing behavior. 

We checked the stability of the obtained solutions by 
performing direct numerical simulations of their dynam- 
ics governed by Eq. , after introducing a perturbation 
of the initial form (1 + 5)4>(x), where <5 was ±0.01 and 
4>{x) was given by one of the distributions shown in Fig. [2] 
(A-C). Performing computations until t = 1000 we did 
not observe any significant change in the behavior of the 
modes In the case (D) , however a decrease of the 
amplitude (S = —0.01) resulted in spreading out of the 
solution, confirming the results reported earlier in [||. 
Increase of the amplitude, achieved by taking <5 = 0.01, 
resulted in a oscillation of the shape of the mode, which 
was neither collapsing nor reaching any stationary state. 



B. The repulsive case. 



Although near the bifurcation point, the shapes of the 
modes collapsing with each other do not allow clear iden- 
tification of n-soliton solutions, the introduced interpre- 
tation is well supported by the analysis of the modes at 
large negative frequencies, illustrated in Fig. [3J Indeed, 
independently on the potential depth, solutions of all 
branches tend to chains of equidistant Townes solitons, 
the number of atoms of each one approaching A* . While 
the lowest branch solution has its counterpart for the ho- 
mogeneous QNLS equation, the coupled solitonic states 
of the higher branches do not exist at V(x) = 0. This 
means that a lattice of arbitrary depth can be treated 
as a singular perturbation allowing soliton binding. Our 
numerical studies allow us to formulate the following con- 
jecture about "quantization" of the number of particles: 

Conjecture 1 - In the limit to — ► -co there exists an 
infinite number of different stationary solutions of the 
QNLS equation with a periodic potential, each of them 
having nN*, where n is an integer, number of particles. 

Recalling that the period of the potential is ir one ob- 
serves form Fig. that the distances between coupled 
Townes solitons are integers of the lattice constant: it is 
2tt and n in panels (b) and (c), respectively. One more 
property of the solutions one can observe in the upper 
panels of Fig. [3J approach of the solutions to a set of 
Townes solitons, as lo — > — oo, is much faster for smaller 
amplitudes of the potential. This is explained by the 
fact, that the respective limit corresponds to the situ- 
ation where the periodic potential can be viewed as a 
(singular) perturbation, and is achieved more rapidly for 
smaller amplitudes of the lattice. 

Summarizing the results illustrated in Figs. |2 and [31 
we concludes that there exists no upper bound on the 
number of particles which can be loaded in lattice. 

Finally, we notice that the whole lower branch in Fig. [21 



In the repulsive case one can easily prove that no lo- 
calized solution of Eq.Q can exist for u> < — Vq. Indeed, 
from the Pohozaev identities for with (7 = 1, which 
read 

J Ul - ~<£ 6 + uxjp" - [V{x] + xV x {x)}^ dx = 
/(^-^ + V<*,^ = 0, 
we deduce 

loN = J (<g + <f> 6 + V{x)ct> 2 ) dx 
> J V(x)<j) 2 dx > -V N 

In the case (|15|) the region where u> < — Vq does not 
cover the semi- infinite gap completely (see Fig^l . In the 
non-covered part of the semi-infinite gap (a finite interval 
of lo existing for any fixed Vo) we performed a search of 
localized solutions using a code based on a specific version 
of the shooting method. The obtained results allows us 
to formulate the following: 

Conjecture 2 - No localized solution of Eq. 0) can exist 
in semi-infinite gap in the repulsive case. 

However we failed to prove the Conjecture 2 rigorously. 

VI. GAP SOLITONS OF THE FIRST GAP. 

The behavior of N vs lu lying in the first gap is shown 
in FigQ] (attractive case, odd modes) and FigOO (repulsive 
case, even modes). As for the semi-infinite gap, one ob- 
serves nonzero minima for the number of particles, which 
are necessary for creation of gap solitons and achieved 
at nonzero detuning toward the gap. Another interest- 
ing feature that one can see on the both figures Fig0| 
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and [3] is that the numbers of particles for the localized 
modes grow infinitely when frequency approaches one of 
the gap edges. This phenomenon is explained by the al- 
gebraic asymptotic of the gap soliton, which corresponds 
to the border of the band: (f> oc l/x 1 / 2 as x — > oo. It 
is worth it pointing out here that, due to the sign of 
the effective mass, which in terms of the Eq. Ill 21) would 
mean aM a Xa > 0, it does not allow existence of small- 
amplitude bright solitons of the type (|13|l in vicinity of 
lower bound of the first gap in attractive case and upper 
bound of the first gap in repulsive case, Vq > 0. The 
following self-consistent recurrence procedure allows one 
to compute the basic term of algebraic asymptotics for 
the gap soliton corresponding to the border of 1-st gap 
as well as higher corrections to it. 

Let us denote by u> a the respective gap boundary, such 
that uj a = u>lt+i h° cr = 1 and u> a — w„ if a = — 1. 
Let (fo be the periodic Bloch function corresponding to 
Lu a (notice that, compared with Sec. IV, here, for the 
sake of convenience, we change the notations of Bloch 
functions bordering a gap), <fo(x) = pa(x + 27r), normal- 
ized by the relation ip^dx = 1. We seek a solution 
of Eq. J3J in the form of the formal asymptotic series 
= 7 (ipo/x 1 / 2 + tpi/x 3 / 2 • • • ) , where 7 is a constant 
to be found. This leads to the recurrent formula for <p n 
(n = 1,2...), the first step of which leads to 



Continuing the described procedure yields the terms of 
the expansion for <j> up to any order and all of them are 
unambiguously defined. Thus, the logarithmic divergence 
of the number of particles N takes place as w — > ii) a , 
which conforms to Figs. [5j Here we notice that the 
similar expansion for the cubic NLS equation gives the 
decay oc x^ 1 , which implies finiteness of the number of 
particles for the solution on the gap boundary. 




FIG. 4: Number of particles of odd modes vs frequency in the 
first gap, for Fo = 3 in attractive case (central box). Mode 
profiles in the points A (u ~ 1-26), B (u> — —0.73), C (lj = 
2.15), and D (u) = 1.33) are shown in the side boxes, ~ 
1.256 is the bifurcation point of the second and third branches. 



dx 2 



+ [u a - V(x)] ifi 



d(fp 
dx 



(16) 



A solution of this equation can be represented in the 
form ipi = ifi(x)+ciifo{x), where a is a constant and (p\ 
satisfies the orthogonality condition (ipi,ipo) = 0. Con- 
sidering the terms of the third order (n — 3), i.e. those 
proportional to x~ 5 ^ 2 , we obtain the equation 

+ [uj a - V(x)\ ip 2 = -0-7 + 3— - -(p (17) 
for which the orthogonality condition yields 

r 2-K n p2tt r 27T 



r 3 r r dip 

/ ifn dx H — / pi dx — 3 / ipo — — dx = 
J a 4 ,/ J dx 



Taking into account the representation of p(x) as well as 
normalization condition for ipo we compute 



12o- f^<p ^dx-3a 
Aj^ipldx 



1/4 



(18) 



Next, representing p>2 = ¥2 + C2P>o, and considering the 
terms corresponding to x~' 7 1 2 we obtain 

d 2 P3 . r ~ „, , } -dip 2 4 4 15 

-^j- + [W a ~ V(x)\ ip 3 = 5— - 57 tp tp! - —ipt (19) 

Now, from the orthogonality condition we obtain the co- 
efficient ci (which is expressed by a bulky formula not 
presented here). 
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FIG. 5: Number of particles of even modes vs frequency in 
the first gap, for Fo = 3 in repulsive case. Mode profiles in 
the points A (w = -0.7105), B (w = -0.73), C (u> = 2.16), 
and D (lu = —0.5) are shown in the side boxes. 

Returning to Fig.^J we observe that all branches of the 
spectrum are located above the critical value N*, which 
suggests the possibility of observing instability and col- 
lapse of the respective solutions. Fig.|SJi shows an exam- 
ple of such dynamics obtained by direct numerical simu- 
lations. 

In both panels presented in Fig. EI the dynamics shows 
three different stages of evolution, the first two being sim- 
ilar for both attractive and repulsive cases. At the first 
stage, one observes growth of the amplitudes and shrink- 
ing of the width of the solutions. This " quasi-collapse" 
behavior is explained by the fact that in both cases the 
dynamics is approximately described by the equation 
(|12fl . which is also an effective QNLS equation, and thus 
corresponds to collapsing solutions. This is not an au- 
thentic collapse, because after the particles are gathered 
mainly in a few (actually two) potential wells, l|12f> is not 
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FIG. 6: Dynamics of gap solitons in attractive (a) and repul- 
sive (b) cases, starting with perturbed initial conditions given 
in Fig.^P and in Fig. 03 multiplied by a factor (1 + 5) with 
S = 0.01. Negative S in the both cases resulted in spreading 
out solutions (not shown). The simulation were carried out 
on the spatial interval (-60,60). 

applied any more and (roughly) the system is similar to a 
dimer model. Oscillatory regime in such a model can be 
proven to originate the dynamics insensitive to the 
sign of the nonlinearity. This is the reason we observe 
qualitatively similar oscillatory behavior of both systems 
during the second stage of evolution. Finally, after sev- 
eral (three in Fig. a) oscillations in the attractive case 
we observed collapse (the third stage). In the repulsive 
case, oscillatory behavior was observed for longer times, 
and the last (third) stage corresponds to decay of the 
peak amplitude and spreading out of the pulse. 

VII. CONCLUSION 

In the present paper we have shown that a periodic 
potential results in qualitative changes in behavior of the 
quintic nonlinear Schrodinger equation. It modifies the 
sufficient condition for the collapse, requiring negative 
energies with larger absolute values, allows the existence 
of infinite families of the localized solutions, and can bind 



Townes solitons (the latter expressed by the conjecture 
about quantization of the number of particles in the limit 
of large negative frequencies). Either by coupling Townes 
soliton or by allowing existence of gap solitons, the pe- 
riodic potential allows storage of an infinite number of 
atoms. It also enforces the stability of Townes solitons 
allowing them to have under-critical number of particles. 

The localized modes were shown to be very dif- 
ferent from their counterparts in the cubic nonlinear 
Schrodinger equation, where the number of particles can 
go to zero and is bounded from above. Localized modes 
in the quintic nonlinear Schrodinger equation on the one 
hand require a minimal nonzero number of particles and 
on the other hand may have arbitrarily large number of 
atoms. 

The results reported in the present work, however, do 
not exhaust all stationary solutions. Speaking about so- 
lutions with of a definite parity we have restricted con- 
sideration only to the families originated by the lowest 
branch (i.e. by the one having the smallest number of 
particles) . Possibility of existence of branches having dif- 
ferent symmetry (say, odd solutions in the semi-infinite 
gap and even solutions in the first gap, in the case of 
attractive interactions) is left to be open question. 
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